Muon spin relaxation and isotropic pairing in superconducting PrOs4St>i2 
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Transverse-field muon spin rotation measurements in the vortex lattice of the heavy-fermion (HF) 
superconductor PrOs4Sbi2 yield a temperature dependence of the magnetic penetration depth A in- 
dicative of an isotropic or nearly isotropic energy gap. This is not seen to date in any other HF super- 
conductor, and is a signature of isotropic pairing symmetry, possibly related to a novel nonmagnetic 
"quadrupolar Kondo" HF mechanism in PrOs4Sbi2. The T = relaxation rate cr s (0) = 0.91(1) /is -1 
yields an estimated magnetic penetration depth A(0) = 3440(20) A, which is considerably shorter 
than in other HF superconductors. 
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In most heavy-fermion (HF) metals and supercon- 
ductors the / ion (Ce, Yb, U) has a magnetic ground 
state. HF behavior has, however, been reported in a 
small number of praseodymium-based alloys and com- 
pounds jl], H H, m which the crystalline electric field 
(CEF) ground state of the non-Kramers Pr 3+ ion could 
be nonmagnetic and degenerate. In this case a charge- 
scattering analog of Kondo spin scattering can give rise to 
the so-called "quadrupolar Kondo effect" j| , an example 
of the two-channel Kondo effect that has been invoked to 
explain non- Fermi- liquid behavior in HF systems ||. 

Superconductivity has recently been discovered in the 
cubic HF compound PrOs4Sbi2 |], @). From thermo- 
dynamic measurements it is found that a large carrier 
effective mass m* ss 50 m e characterizes both the normal 
and superconducting states, and the transition tempera- 
ture T c = 1.85 K is relatively high for a HF superconduc- 
tor. Although a conventional spin-based HF mechanism 
has not been completely ruled out, thermodynamic prop- 
erties of PrOs4Sbi2 suggest a Pr 3+ nonmagnetic doublet 
T$ CEF ground state, so that the quadrupolar Kondo ef- 
fect is a candidate mechanism for the HF behavior. The 
symmetry of the superconducting pairing in such a sys- 
tem is a fundamental question. 

Transverse-field muon spin rotation (TF-//SR) has 
proved an effective probe of the internal magnetic field 
distribution n(B) in the vortex state of conventional 
and unconventional type-II superconductors || |[ [To). 
In the ^SR technique @ [l|] spin-polarized positive 
muons are implanted into the sample, where the 

muon spins precess in their local fields. In general 
the muon spin relaxation rate is related to the rms 
width [(AS) 2 ] 1 / 2 of n(B). In turn [(AS) 2 ] 1 / 2 is inversely 



proportional to the square of the magnetic penetration 
depth A, which is related to the density n s of supercon- 
ducting carriers and m* by the London equation 

1/A 2 = 4vrn s e 2 /m*c 2 . (1) 

The temperature dependence of A at low temperatures 
is therefore sensitive to the lowest-lying superconduct- 
ing excitations, the thermal population of which reduces 
n s with increasing temperature. In the presence of an 
isotropic or nearly isotropic energy gap A(T) — A(0) varies 
exponentially with temperature, whereas nodes in the 
gap function characteristic of non-s-wave pairing lead to 
power-law dependences || ^, [lO] . In HF superconductors 
A is typically very long (> lCr A Q), in which case the 
muon relaxation rate can be dominated by other sources 
of magnetism in the sample. It is then difficult to extract 
A(T) from TF-/L/.SR experiments, but the technique has 
been successfully applied in a number of systems Q. 

This Letter reports /iSR experiments in the super- 
conducting state of PrOs4Sbi2- The penetration depth 
derived from the vortex-lattice field distribution width 
exhibits the temperature dependence characteristic of 
isotropic s-wave pairing. Isotropic p-wave pairing 
which is indistinguishable from s-wave pairing by ther- 
modynamic or electrodynamic measurements, is also pos- 
sible. To our knowledge this is the only example to date 
of an isotropic gap in a HF superconductor. It suggests 
the possibility of (a) marked differences in superconduct- 
ing properties between HF materials with magnetic and 
nonmagnetic /-ion ground states, and (b) a relation be- 
tween the pairing symmetry and the mechanism (spin or 
quadrupole Kondo effect) for the HF normal state from 
which the superconductivity evolves. 
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/iSR measurements were carried out in the dilution re- 
frigerator at the M15 beam line, TRIUMF, Vancouver, 
Canada, on an unaligned powder sample of PrOs4St>i2 
prepared by growth from an antimony flux as described 
previously || |7j. Gallium arsenide, which rapidly depo- 
larizes muons in transverse field, was placed around the 
sample to minimize any background signal due to muons 
that stop in the sample holder. TF-/1SR and zero-field 
/iSR (ZF-/1SR) data were obtained for temperatures be- 
tween 0.05 and 2.1 K. 

In a number of HF superconductors local fields due to 
spontaneous static magnetism below T c are inferred from 
ZF-/1SR data H 0]. Such magnetism can originate from 
a spin-wave instability or local-moment ordering, as ob- 
served in CeCu2.2Si2 H where magnetic ordering and 
superconductivity compete for the ground state. Alter- 
natively, static magnetism below T c may be due to a su- 
perconducting order parameter that breaks time-reversal 
symmetry, as observed in, e.g., Th-doped UBei3 M, [l5| . 
Muon local fields due to static magnetism could mas- 
querade as vortex lattice inhomogeneity, and falsify the 
interpretation of the TF-/iSR results given above. 

We have therefore carried out ZF-/iSR experiments 
to determine whether such static magnetism exists in 
PrOs4Sbi2. The zero- field data are well fit by a "damped 
Kubo-Toyabe" function, i.e., the product of a damping 
exponential and the Kubo-Toyabe function Q expected 
from nuclear dipolar fields. Figure [j] gives the temper- 
ature dependence of the corresponding relaxation rates 
W and Akt- The exponential damping rate W(T) is 
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FIG. 1: Damped Kubo-Toyabe ZF-^tSR relaxation in 
PrOs4Sbi2. Circles: Kubo-Toyabe relaxation rate Akt(T). 
Triangles: exponential damping rate W(T). 

already appreciable in the normal state, increases only 
slightly more than experimental uncertainty below T c , 
and is (negatively) correlated with Akt (I 1 ) in the fit- 
ting process. Thus there is no statistically significant 
evidence for static magnetism below T c ; the present data 
place an upper limit of ~ 50 [iT on any static field below 
T c . We argue below that the relaxation rate changes of 
Fig. |] are too small to affect the considerably increased 



TF-//SR rates in the superconducting state. 

Representative TF-^SR muon-spin precession signals 
at an applied field of 20 mT are shown in Fig. || in the 
normal and superconducting states. A small (~10%) 
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FIG. 2: TF-/1SR spin precession signals in PrOs4Sbi2 (T c — 
1.82 K), applied field 20 mT. (a) Normal state (T = 2.1 K). 
(b) Superconducting state (T — 0.1 K). The weak nonre- 
laxing signal in (b) is due to muons that do not stop in the 
sample. 

nonrclaxing background signal is visible at long times in 
Fig. ||(b). Normal-state data for temperatures just above 
T c are best fit by an exponential decay function, with a 
rate constant W n — 0.17(1) /is -1 at 20 mT, rather than 
the Gaussian commonly found when the broadening is 
due to nuclear dipolar fields ]l0| . This suggests that the 
broadening is due to a different mechanism, most likely 
inhomogeneity in the Pr-moment susceptibility. Consis- 
tent with this view, in the normal state W n grows linearly 
with applied field (data not shown). 

The internal field distribution in the vortex state is the 
convolution of the field distributions due to the vortex 
lattice and to the electronic and nuclear moments of the 
host material; by the convolution theorem the muon-spin 
precession signal is the product of the Fourier transforms 
of these distributions. The superconducting-state data 
have been fit to a product of exponential and Gaussian 
functions, where the latter is taken as an approximation 
to the relaxation function due to the vortex-lattice field 
distribution. The exponential rate W n was held fixed at 
the normal-state value, although at low temperatures the 
values of the superconducting-state Gaussian rate a s are 
not greatly changed if W n is allowed to vary for best fit. 

The fit is statistically satisfactory, as can be seen qual- 
itatively in Fig. 0(b). This indicates that n(B) is essen- 
tially Gaussian, due perhaps to disorder in the vortex 
lattice; the non-Gaussian relaxation expected from the 
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field distribution in an ideal vortex lattice M, |17|] was 
not observed. The Gaussian relaxation rate a s is then 
given by 7m [(AB) 2 ] 1 / 2 , where 7^ = 8.516 x 10 8 s" 1 T" 1 
is the fj, + gyromagnetic ratio. 

The temperature dependence of a s is shown in Fig. || 
for an applied field of 20 mT. Salient features of these 

PrOs 4 Sb 12 ji Q H= 20 mT 
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FIG. 3: Temperature dependence of vortex-state n + relax- 
ation rate cr s (T) in the superconducting state of PrOs4Sbi2. 
Curve: a s (T) = a s (0)[l - (T/T c ) v ], a a {0) = 0.91(1) /is' 1 , 
T c = 1.83(2) K, y — 3.6(2). Inset: low-temperature penetra- 
tion depth A(T) derived from a s (T). Curve: A(T) = A(0)[1 + 
(7rA/2T) 1/2 exp(-A/T)], A(0) = 3440(20) A, A/T c = 2.1(2), 
from fit to data for T < 0.8 K. 

data are that cr s (T) varies only slowly at low tempera- 
tures, and that the T = value a s (0) = 0.91(1) /is -1 . 
This is large for a HF superconductor ^ ^) and is much 
larger than changes in the ZF-/iSR rates (Fig. [l]) , so that 
a s is dominated by the vortex-state field distribution. 
We obtain A from a s and the expression 



(Ai?) 2 = 0.00371 ^A -4 



(2) 



appropriate to an isotropic extreme type-II superconduc- 
tor ]l7| , where $0 is the flux quantum. The data are 
shown in the inset to Fig. [| The BCS low-temperature 
expression A(T) = A(0)[1 + (ttA^T) 1 / 2 exp(-A/T)] JTJ 
was fit to data for T < 0.8 K < 0.4 T c ; the curve in the 
inset to Fig. || shows this fit and its extension (dashed) 
up to 1.2 K. The fit value A(0) = 3440(20) A is short for 
a HF superconductor. The fit value of the ratio A/T c = 
2.1(2) is somewhat larger than the BCS value 1.76, 
suggesting strong coupling. Over the entire temper- 
ature range cr s (T) is consistent with the phenomeno- 
logical "two-fluid" temperature dependence 1/A 2 (T) oc 
1 — (T/T c ) 4 , although the data are slightly better fit with 
an exponent of 3.6(2) (curve). All these properties in- 
dicate that the gap is isotropic or nearly so ||, ||, [uj. 
Proportionality between (AB) 2 and A~ 4 should survive 
vortex lattice disorder as long as the distance between 
vortices is much smaller than A (H ^> H c i); disorder in- 
creases the numerical coefficient in Eq. (||) Jl7| but should 
not affect the temperature dependence of (AB) 2 . 



Early studies of the penetration depth in cuprate su- 
perconductors using unaligned powders also observed lit- 
tle temperature dependence of A at low temperatures, 
and concluded that the pairing was s-wave |jjo[ . Later 
measurements in high-quality aligned crystals fpfi|] re- 
vealed the linear temperature dependence characteristic 
of (i-wave pairing. The early results appear to have been 
due to a combination of circumstances: strong anisotropy 
in A, oxygen inhomogeneity, and sensitivity to hole dop- 
ing due to the low density n s of superconducting carriers. 
None of these factors would seem to affect the current 
measurements. The crystal structure of PrOs4Sbi2 is 
cubic, so that the penetration depth is isotropic. With 
to* w 50m. e we find n s 10 22 carriers/cm 3 from Eq. (Q), 
so that PrOs4Sb 12 is a good metal and the carrier con- 
centration is insensitive to chemical inhomogeneity. Thus 
it seems unlikely that the observed temperature depen- 
dence of A(T) is due to extrinsic effects. 

The resistive mean free path I can be calculated from 
the residual resistivity p(0) ~ 5 /xf2-cm R using n s 
as an estimate of the normal-state carrier concentration, 
and can be compared with the superconducting coher- 
ence length £ @. We find £/£ w 3, so that PrOs 4 Sbi 2 
is a rather clean superconductor [and Eq. (Q), derived 
in the clean limit, is valid]; this is of importance in the 
analysis of thermodynamic properties. 

The field dependence of a s at T = 0.1 K in PrOs4Sbi2 
is shown in Fig. 0. In these experiments the field was 
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FIG. 4: Dependence of TF-/iSR relaxation rate a a (H) on 
applied field in PrOs 4 Sbi 2 , T = 0.1 K. 

set with the temperature above T c , and the sample was 
cooled into the superconducting state at constant field; 
otherwise trapped flux broadens the field distribution 
and causes a spurious increase of the /i + relaxation rate. 
The field dependence of cr s is nonlinear, decreasing at 
low fields and then saturating above ~ 100 mT. The 
origin of this nonlinearity is not understood, but might 
be due to disorder in the vortex lattice |l0|, or to a 
change of vortex-lattice symmetry, e.g., from a triangu- 
lar to a square lattice. Measurements at fields higher 
than ~ 200 mT will be difficult because of the growing 
importance of the normal-state paramagnetic relaxation 
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rate. 

Figure || is an Uemura plot (if], ||(J of superconduct- 
ing transition temperature T c vs zero-temperature muon 
relaxation rate cr s (0) in a number of HF superconduc- 
tors: UBei3 f2lfl , CeCu2Si2 (derived from the penetration 
depth estimate of Ref. §), UPd 2 Al 3 §§, CeCoIn 5 and 
Celrln 5 , [||, UPt 3 §1, and PrOs 4 Sbi 2 (this Letter). It 
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FIG. 5: Uemura plot of superconducting transition tem per 
ature T c versus muon relaxation rate a s (0) for UBei3 
CeCu 2 Si2 §, UPd 2 Al 3 @, CeTIn 5 , T = Co and Ir ||| 
UPt 3 §J, and PrOs 4 Sbi2~(tliis Letter). The "Uemura line' 
is followed by a number of unconventional superconductors. 



can be seen that a s {0) is larger in PrOs4Sbi2 than in 
other HF superconductors. The values of cr s (0) shown in 
Fig. ^ arc in rough inverse proportion to the HF effective 
masses of these systems, as expected from Eq. (|l|) and 
cr s (0) oc 1/A 2 (0), so that the position of PrOs4Sb 12 in 
this plot is not surprising. 

The isotropy of the superconducting energy gap in- 
dicated by the temperature dependence of a s (Fig. ||) 
strongly suggests that superconductivity in PrOs4Sbi 2 is 
in some sense "conventional." In this regard PrOs4Sbi 2 
differs markedly from other HF superconductors stud- 
ied to date, a difference that may possibly be related 
to a nonmagnetic or quadrupolar Kondo state in this 
compound. Better understanding of the quadrupolar 
Kondo lattice will be needed to elucidate the questions 
of whether PrOs 4 Sb 12 is indeed such a system and, if 
so, whether this property is related to the conventional 
superconducting behavior. 
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